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ABSTRACT 

A single magnetic monopole in pure SU(2) gauge theory is simulated on the lattice 
and its mass is computed in the full quantum theory. The results are relevant for 
our proposed realization of the dual superconductor hypothesis of confinement. 

1. Introduction 

In Ref. 1 we presented a realization of the dual superconductor hypothesis of con- 
finement in S U(2) gauge theory. The monopoles were of the 't Hooft-Polyakov (HP) 
form, in particular the Bogomol'nyi-Prasad-Sommerfield (BPS) solution. Assuming 
that the configurations considered dominate the string tension a we obtained the 
estimate v^/^ms = 2-3. 

In our model we made two assumptions. The first one was that quantum fluc- 
tuations screen the (chromo-)electric field of the classical BPS monopole, turning it 
into what we call an HP-like monopole. The second assumption concerned the mass 
of the monopole. The mass of the classical BPS monopole of scale ji is M = Airfi/g 2 . 
We conjectured that the monopole mass in the full quantum theory behaves as 

M = 2 ^ , C(gUKR^)) • (1) 

Here g R is the running coupling in the i?-scheme defined in terms of the quark- 
antiquark potential, with Ar/Aj^ = 1.048. The function C was supposed to be 
slowly increasing as for the HP monopole, 1 < C ^ 2. 

In the present study, described in detail in Ref. 2, we investigate the second as- 
sumption numerically. We put a single monopole in a box, by imposing appropriate 
monopole boundary conditions, and use lattice Monte Carlo methods to include 
quantum fluctuations. The simulation data are used to determine the function 
C(g R (A R / fi)) which allows us to check the validity of the mass formula (1). 

2. The magnetic monopole 
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The classical magnetic monopole of (arbitrary) scale fi in euclidean pure SU(2) 
gauge theory is given by the gauge field configuration 



A a k {x, t) = e aklXl ^ , A a 4 (x, t) = S ak x k ^^ 



with 

cosh fir fir 

H(fj,r) = fj,r-— 1, K(fj,r) = -— , (3) 

smn fir sinn fir 

as in the BPS limit of the HP monopole. 

The idea is to induce the monopole in a finite cubic spatial box (times a periodic 
time direction of length T) by fixing the fields in the boundary at a value suggested 
by the asymptotic behaviour of the monopole field. We take 'HP-like' boundary 
conditions with parameter fi for our dynamical simulations, i.e. K ~ 0, H/r ~ fi . 
An analysis of classical monopole energies shows that this is compatible with both 
an HP monopole of scale fi = fi and a BPS monopole of scale /i e ff = fi + l/i? e ff- 
Another important result of the classical analysis is that the accessible range of fi 
values is restricted to < fio < ir/T. This is due to a symmetry implying that fi 
and fi' boundary conditions are equivalent if fi' = fi + 2irn/T } for some integer n, 
and to the occurrence of monopoles of opposite electric charge. 



3. The monopole mass 



The monopole mass can be written as the sum of a contribution from inside the 
box, 

measured in the simulation, and a correction term for the outside region. Here 
Z mon ((3;afi ) = jDUexp[-S a (U;(3)} = J DU exp[-(3 £ (1 - ^tr □(£/))] (5) 

plaquettes 

□ 

is the partition function subject to monopole boundary conditions and Z VSLC ((3) is 
the analogous definition for vacuum boundary conditions, i.e. = 0. The gauge 
coupling constant g is embodied in fi = 4/g 2 and fi — > oo is the classical limit. 

Eq. (4) can be written in a form accessible to Monte Carlo computation by 
differentiating it with respect to fi and subsequently integrating again, 

M m (ft) = f^E(fi) = - (^)vac) • (6) 

Jo fi Jo fi i 

In order to compute this numerically the integral is replaced by a sum. At each 
value of fi in this summation two simulations are needed, to compute (S) with both 
monopole and vacuum boundary conditions. High statistics is required to compute 
the difference of these two large numbers accurately. 



Subsequently, C- m (g 2 R (hnl ji)) is extracted from the M- m data as follows. First 
Ar/ fi is calculated using Monte Carlo data for ay/a and taking a value for y/a/ Ar. 
Next q r {ArI ji) is calculated using its two-loop /^-function. Inserting this into Eq. (1) 
leads to C- m (g R (AR/ //)). This procedure is carried out for both ji = /i and ji = /i e ff, 
corresponding to HP-like and BPS-like behaviour of the monopole, respectively. 

Fig. 1 shows C for a simulation at an 8 4 lattice, with fi = ir/8 = 0.39. The 
correction due to the mass contribution from the outside region is included. For 
g R mO we know the monopole is BPS-like, so the HP analysis (upper three curves) 
is misleading there. At larger couplings, the monopole may or may not become 
HP-like. If it does, the C values will lie in the region indicated by the upper set 
of curves. This means that C increases from C = 1 at weak coupling to C ~ 1.6 
a ^ 9r ~ 8 or C ~ 2.0 at g R 10, depending on the value of \J~a / Ar. This is in 
good agreement with our assumptions. If, however, the monopole remains BPS-like 
at large coupling, our first assumption does not apply. Nevertheless, even in that 
scenario (lower set of curves) C remains almost constant, C ~ 1 up to g R 6. 
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Fig. 1. C as a function of (/|j(A/j//i). The different sets of points are for \fa / Ar = 1.7 (O), 2.0 
(□) and 2.2 (A), for both jl = [1q (upper set of curves) and jl = /i e ff (lower set). The horizontal 
lines denote the classical limit (g 2 = 0). The solid curves come from large-/? fits of the mass data. 
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